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Abstract 

It is well known that in string compactifications on toric Calabi-Yau manifolds one 
can introduce refined BPS invariants that carry information not only about the charge 
of the BPS state but also about the spin content. In this paper we study how these 
invariants behave under wall crossing. In particular, by applying a refined wall crossing 
formula, we obtain the refined BPS degeneracies for the conifold in different chambers. 
The result can be interpreted in terms of a new statistical model that counts 'refined' 
pyramid partitions; the model provides a combinatorial realization of wall crossing and 
clarifies the relation between refined pyramid partitions and the refined topological vertex. 
We also compare the wall crossing behavior of the refined BPS invariants with that of the 
motivic Donaldson-Thomas invariants introduced by Kontsevich-Soibelman. In particular, 
we argue that, in the context of BPS state counting, the three adjectives in the title of 
this paper are essentially synonymous. 
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1. Introduction 

This paper is devoted to ttie study of ttie space of BPS states, H-bps^ in type II 
string compactifications on Calabi-Yau 3-folds. In general, such compactifications lead to 
effective Af = 2 theories in four dimensions, and, by definition, Hbps is the subspace 
of the Hilbert space of an effective four-dimensional theory that consists of one-particle 
states transforming in small representations of the d = 4, Af = 2 supersymmetry algebra. 
The space Hbps encodes much interesting information about the Calabi-Yau space X as 
well as about the physics of the four-dimensional M = 2 theory, in particular providing 
connections to black hole physics and topological strings [Q. 

Let us summarize some of the basic properties of the space H-bps- First, it is graded 
by charge sectors, 

Hbps = 0Hbps(7), (1-1) 

where F denotes the charge lattice. For example, in type IIA string theoryS on a Calabi- 
Yau 3- fold X the BPS states in question are bound states of Dp-branes with even values 
of p, so that F = Z) and 

-f = ch{S)/AiX) = pO + p + Q + 

DQ DA D2 DO 

is the charge vector of a D6/D4/D2/D0 bound state ( equivalent ly, the Mukai vector of the 
corresponding coherent sheaf S). Roughly speaking, in this case Ti-BPsil) — H* 
where AA{'y) denotes the moduli space of branes of charge 7. In addition, Hbps a, rep- 
resentation of the rotation group Spin{3) in four space-time dimensions. This gives Hbps 
an extra grading, which eventually leads to the refinement of BPS invariants considered 
below. Thus, altogether, Hbps comes equipped with a F © Z-grading. 

The space of BPS states, Hbps, depends moreover on the asymptotic boundary con- 
ditions in four space-time dimensions. Much of its interest actually comes from the de- 
pendence on this extra data, which includes the moduli of the Calabi-Yau 3- fold X f^. 

^ As is well known, type IIA string theory on a Calabi-Yau space X is dual to type IIB string 
theory on a mirror Calabi-Yau space X. In what follows, we pick a duality frame corresponding 
to type IIA theory. 
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Instead of working directly with Hbps, it is often convenient to consider a simpler object, 
the index of BPS states! 

n{r,u):=Ttni-y;u){-lf (1-3) 

that 'counts' BPS states of given charge 7 and is invariant under complex structure de- 
formations of X. However, as the notation indicates, the index 0(7; u) still depends on 
the asymptotic value of the complexified Kahler moduli, u = B + iJ. It is a piecewise 
constant function of u that can jump across walls of marginal stability, where the phases 
of the central charges of the constituents of a bound state align. 

Our main focus in this paper will be the refined BPS index, defined a^ 

(7; u; y) = J2{-y)-nf{r, u) := T^n^^.u) {-yf' , (1.4) 

n 

where J3 is a generator of the rotation group Spin{3). In a simplified situation where 
T^BPsin'iU) admits a description as the cohomology of the brane moduli space A^(7;w), 
the BPS index 0(7; u) and its refinement 0''*^^(7; w; y) correspond, respectively, to the Euler 
characteristic and to the Poincare polynomial of A^(7; u). In particular, from the definition 
( p. .41) it is clear that at y = 1 we have 0'''^^(7; u; 1) = 0(7; u). 

While the refinement 0'^'^^(7; u; y) captures useful information about the spin content 
of BPS states, it is generically not invariant under complex structure deformations of X. 
However, if X has no complex structure deformations, then the refined BPS invariants 
0'''^^(7; u] y) as well as the space Hspsil', u) itself are expected to be interesting invariants 
of X . Thus, in certain examples the refined BPS invariants can be related to equivariant 
instanton counting or to categorification of quantum group invariants 0. In fact, 
in the early days of the refined BPS invariants, the only practical way to compute them 
was by using one of these relations. The situation improved significantly with the advent 
of the refined topological vertex which reduced the computation of refined BPS invari- 
ants for an arbitrary non-compact toric Calabi-Yau 3-fold X to a systematic combinatorial 

^ More precisely, this index is the second helicity supertrace. Following in our definition of 
Tin PS we tacitly factored out the contribution of a universal half-hypermultiplet associated with 
the position in R^, allowing us to write the index of BPS states in the simple form given here. 
In TLbps, a hypermultiplet counts as a state of spin zero (hence = 1) and a vector multiplet as 
spin 1/2 (hence Q = —2). 

Again we factor out a contribution y^^(l — y)^ from a universal half-hypermultiplet. 
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algorithm based on the counting of 3D partitions. Note that all toric Calabi-Yau mani- 
folds are automatically rigid and provide an excellent laboratory for studying refined BPS 
invariants. They will also be our main examples in the present paper. 

In the context of local toric Calabi-Yau manifolds, a natural object to consider is a 
generating function 

Z{q,Q;u):= {-qTQ^ ^hp,n;u) (1.5) 

that 'counts' BPS states of DO and D2 branes bound to a single D6 brane. Here, 'yp,n is a 
shorthand notation for the charge 7 = (1,0, —(3, n) of a D6/D2/D0 system with n units of 
the DO-brane charge and the D2-brane charge corresponding to a curve in homology class 
P G H2{X; 2). In one of the chambers, the D6/D2/D0 partition function (|1.5| ) is the usual 
generating function of Donaldson-Thomas/Gopakumar-Vafa invariants PJ^I- Similarly, the 
refinement of (|1.5| ), 

(q, Q, y- u) := i-^yQ^ ^'""'hp^n, u- y) , (1.6) 

/36H2(^;2) 

carries information not only about the charges of the D6/D2/D0 bound states but also 
about the spin content. For u in the DT region of the Kahler moduli space, it reduces 
to the generating function of the usual refined BPS invariants computed by the refined 
topological vertex. One of our goals will be to study the refined partition function ( pT^ ) 
in other chambers and, more generally, to understand how the refined BPS invariants 
0'^'^^(7; u] y) change across walls of marginal stability. In particular, in the Szendroi region 
we obtain a refinement of the non-commutative Donaldson-Thomas partition function |]^ . 
For example, for the resolved conifold X = Cpi(— 1) © Opi(— 1), it looks like 

00 

^ncdt(9i, 92, Q) = M{q^, q^f J] " Q^i~'^ ^2~'^) (l " Q~' Q^2~~') . (1-7) 

where instead of q and y we use the variables gi and q2, standard in the literature on the 
refined BPS invariants, cf. 

q 

Qi=Qy, q2 = -- (1.8) 
y 

For toric Calabi-Yau 3-folds in the Szendroi region of moduli space, the generating 



functions of BPS states can be computed via statistical crystal melting models |8,10,11 



that are seemingly distinct from the topological vertex formalism [12,13]. Such models 
can be derived from quivers and brane tilings pi]. In the case of the conifold, the sta- 
tistical models take the form of "pyramid partitions," and it was shown that they 
can be generalized to describe many chambers in the resolved conifold moduli space out- 
side the Donaldson-Thomas region. In analogy with the refinement of the topological 
vertex, there exists a refinement of pyramid partitions that computes (for example) the 
generating function ( |1.7| ) . A combinatorial shuffling operation on pyramid partitions corre- 
sponds to refined wall crossing between chambers of moduli space, and, as one approaches 
the Donaldson-Thomas region, the refined pyramid partitions actually resolve into refined 
topological vertices. We expect this behavior to be quite general for toric Calabi-Yau 
manifolds. 

We observe that the wall crossing behavior of the refined BPS invariants is very close 
to that of the motivic Donaldson-Thomas invariants defined by Kontsevich and Soibel- 
ITBI. In fact, the motivic Donaldson-Thomas invariants can also be viewed as a 



man 



"refinement" of the numerical Donaldson-Thomas invariants that depends on the extra 
variable L (the motive of the affine line). Equivalently, this variable can be interpreted as 
a "quantum" deformation parameteiS q in quantization of the complex torus Tr = F^^C*, 
a fact that was extensively used in We claim that this is not an accident and the 

refined BPS invariants are actually the same as the motivic BPS invariants of Kontsevich 
and Soibelman, provided we identify L (resp. q) with the extra variable y"^ that appears 
in the definition (|1.4|) of the refined BPS invariantsi: 

Refined Motivic Quantum 

y „ „ -q1/2 ^1-9) 

Further evidence for this identification comes from the connection with homological in- 



variants of knots and 3-manifolds that will be discussed elsewhere [0. Note that the 
semi-classical limit corresponds to q^^"^ —1 (resp. y — * 1). 

As is well know, in string compactification on a Calabi-Yau 3-fold, the variable y can 
be interpreted as a graviphoton background in four space-time dimensions. Therefore, we 
propose the following 



^ Not to be confused with the formal variable in the generating function Z[q, Q; u). 

^ We are intentionally a little imprecise here. We intend to say that the deformation by y on 
the 'refined' side corresponds to the deformation by —q^^'^ on the 'motivic' side. While in general 
there may be a non-trivial map between these two deformations, as we explain below in section 3 
these deformations actually agree to the leading order. 
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Conjecture: In string theory on a Calabi-Yau space X, the q- deformation of ]7q/ corre- 
sponds to turning on a graviphoton background (a.k.a. Q-background) on R^. In particular, 
we have 



0-f(7;«;y) = 0-°n7;«) 



(1.10) 



with the appropriate identification of variables 

Although in practice it is easier to work with numerical invariants such as 0(7; u) or 
u; y), we should emphasize that the ultimate goal is to "categorify" the wall crossing 
formulae for the numerical BPS invariants and to explore the properties of the space Ti-BPS 



itself. We believe this should lead to a rich mathematical structure. (See [ p!8| , p!9| for earlier 
work where the homological algebra of BPS states played an important role.) In fact, 
it seems that the homological algebra of Hbps is unavoidable if one tries to study the 
refined BPS invariants. For example, there can be walls — we call them "invisible walls" 
— where the ordinary index 0(7; m) doesn't change, but Q^'^^ {■y;u;y) jumps. The basic 
mechanism for how this happens can be understood in a simple situation where H-bps 
admits a description as cohomology of the brane moduli space A4. Then, as u crosses an 
invisible wall, A4 can develop a singularity and undergo a topology-changing transition, 
so that the Poincare polynomial of Ad changes while the Euler characteristic does not. 



These invisible walls have been observed in 19 in certain instances of BPS state 



counting and will be discussed further in [17]. One of their most interesting properties is 



that the best way to describe the change in the spectrum of refined BPS invariants across 
them is by observing that in Kahler moduli space these walls are located where two BPS 
states in short multiplets can combine into a long multiplet. Therefore, the states that 
disappear from the spectrum are trivial in cohomology of the BRST differential 

Q-.HIps^H'^bVs^ (1-11) 

where n = 2J3 denotes the 2-grading of Hbps by spin. (Note that Q changes spin by ^.) 
At the level of the generating function (|1.6| ), it means that across each invisible wall the 
change in Z^'^^{q, Q, y; u) includes an elementary factor {1 — y). 

Organization of the Paper 

In section 2, following (see also ||T5| , pO| , pT| ) , we briefly review some relevant facts 



about the wall crossing behavior of the BPS index 0(7; u) and how it can be used to com- 
pute the D6/D2/D0 partition function in concrete examples, such as the resolved conifold. 
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We then generalize this discussion to refined BPS invariants and compute Z^^^ {q,Q,y;u) 
for the conifold in different chambers, the refined non-commutative Donaldson-Thomas 
partition function ( |1 . 7| ) being a special case. In section 3, we compare the wall crossing 
behavior of the refined BPS invariants with that of the motivic Donaldson-Thomas invari- 



ants introduced by Kontsevich-Soibelman ||T^. In particular, one way to see the proposed 
identification ( |1.1C1|) is to deduce the primitive (or semi-primitive) wall crossing formula 
for the refined BPS invariants from the motivic wall crossing formula of Kontsevich and 
Soibelman. As we shall see in section 3, in this derivation, one is naturally led to the iden- 
tification of the "deformation" parameters in eq. ( |1.9[ ). Finally, in section 4, we interpret 
the results of section 2 in terms of the new statistical model that counts refined pyramid 
partitions. In both refined and unrefined cases, we explain how a shuffling operation on 
pyramid partitions provides a combinatorial realization of wall crossing, and we clarify the 
relation between pyramid partitions and the refined and unrefined topological vertices. 



2. Refined Wall Crossing 

The main goal of this section is to study the wall crossing behavior of the refined BPS 
invariants 0'^'^^(7; m; y) and, in particular, to compute the refined D6/D2/D0 partition 
functions Z^'^^[q^Q^y]u) in different chambers. 

Walls and Chambers 

We are interested in walls of marginal stability for decays 

7 ^ 71 + 72 ■ 

In the Kahler moduli space, such decays take place at the points u where the central 
charges Z{^i]u) and 2(72;^) of the constituents align; the corresponding walls will be 
denoted as W(7i,72): 

>V(7i, 72) = { w I 2(71; u) = AZ(72; u) for some A G IR+} . (2.1) 

Note, in particular, that all walls W(A^i7i, A^272) coincide. 



>V(7i, 72) = >V(Ari7i, N2-f2) , A^i, Ar2 G 2+ 
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The bound state with charge 7 = 71 + 72 is stable on the side of the waU W(7i, 72) where 



(71, 72) ImZ(7i; w)Z(72; u) > . (2.2) 

This condition, as weU as the position of the waU ( ^.11 ), has an elegant interpretation in 
supergravity, if all the charges 7^ are large [H . For example, in the case of a 2-center bound 
state, the separation of the two constituents in is determined by their charges, 

D 1 / \ 1^1 + ^2|cx) 

-K12 = -(7I572): 



2^' ' ' 'lm{ZiZ2] 



IS 



where 2(7) is the central charge function. In particular, we see that the condition ( 2.2 
necessary for the distance R12 to be positive, and R12 diverges when u approaches a wall 
of marginal stability VV(7i, 72). For future reference, we also note that the 2-center bound 
state carries angular momentum [|[: 

Ji2 = ^ (1(71,72)1-1). (2.3) 

In the context of counting BPS states of DO and D2 branes bound to a single D6 
brane on a toric Calabi-Yau manifold X, one encounters walls of marginal stability with 
7i = (1, 0, — m', n') and a "halo" of particles or "fragments" of charge 72 = (0, 0, —rrih, Uh)- 
Following [0,0, we denote such walls as W^'". For example, for the resolved conifold 



X = Opi{—l) © (9pi(— 1), the non- vanishing Gopakumar-Vafa invariants ^2 

0(7 = (0,0,±l,n)) = 1 nGZ 

(2.4) 

0(7 = (0,0,0, n)) = -2 n^O, 
imply that the only walls are 

: D2/D0 fragments 

: m/DO fragments (2.5) 
W° : -DO fragments . 

This leads to the following picture of walls and chambers (in the one-dimensional space 
parameterized hj (f) 0: 



"PV'' w" 



o 
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DT 



core 



-7 
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Fig. 1: The picture of walls and chambers for the resolved conifold proposed in Q. 



(2.6) 



We follow the conventions of Qj^? and denote the resulting chambers for the resolved 
conifold as 

Wall Crossing Formulae and Refinements 

If both 7i and 72 are primitive, the states lost from Ti.BPsili'^) are [|2| 

A^BP5 = (^12) ® ^(71; Ums) ® ^(72; w^s) , (2.7) 



where J12 is the spin of the bound state (p.3|). The corresponding change of the BPS index 
( pT3| ) is given by 

This is the so-called primitive wall crossing formula. It has a very explicit, but lengthier, 
generalization in the case when one of the charges is not primitive. It is then convenient 
to arrange the answer in a generating function that allows one to describe the change of 
the spectrum of BPS states as a sum over Fock spaces [0: 

0x^^A7^|^^^,+Ar,^, =7^ (71) (g)-F(x'=( J^^fc^J ® 7^(^72)) , (2.9) 

N2 k 

SO that 

J^(7i) + An(7i + iV72)x^ = 0(71) (-l)'=<Ti'^2>^fc)fc|(7i,72>|n(fc72) _ ^2.10) 

N>1 k>l 



This is the semi-primitive waU crossing formula. Note that • •)U=o gives the primitive 
wall crossing formula ( |2.8| ). 

In general, if neither 71 nor 72 is primitive, the change of the BPS index across a wall 
of marginal stability is given by the Kontsevich-Soibelman wall crossing formula [ p!B| , pB| 
that will be discussed in section 3 below. We note in passing that if one is interested only in 
counting D6/D2/D0 bound states on toric Calabi-Yau manifolds, then the primitive ( |2.8| ) 
and semi-primitive (|2.10| ) wall crossing formulae and their refinements (which we discuss 



momentarily) are sufficient to determine the D6/D2/D0 partition functions. 

Now, let us generalize this discussion to the refined BPS invariants that, along with the 
charges of BPS states, also carry information about the spin content. Again, starting with 
the simplest case when both charges 71 and 72 are primitive, the change of Q,^'^^{'y,u;y) 
defined in ( |1.4|) is given by 

An^^\r.u;y) = ^ ^\ . ; , 0=^^f(7i;«;2/)0-f(72;t.;y) 

{-y)-{-y) ^ 

= + l + (-y)-<7i,72)+3 ^ _ _ _ ^ (-y)<'/i'7^)-l)0-f(7i;«;y)0-f(72;«;y) . 

(2.11) 



This is the refined version of the primitive wall crossing formula (2^), whose explicit 
form appeared in Indeed, note that setting y = 1 gives the primitive wall crossing 
formula 



For our purposes, we need a refinement of the semi-primitive wall crossing for- 
mula ( |2.1UD . It takes the form 



f^"'(7i)+ 5^ AO=^^f(7i+iV72) 

iV>l 

fc|(7i,72>| (2.12) 

fc>l j=l n 



where are the coefficients of the refined BPS index (|1.4| ). This refinement of ( |2.10| ) 



describes the change of the spectrum ( p.9[ ) while keeping track of the spin of BPS states, 
and satisfies two obvious requirements: 
i) at y = 1 it reduces to the ordinary semi- primitive wall crossing formula ( |2.1C1| ); 
ii) - ■ •)U=o gives the refined primitive wall crossing formula ( p.ll|) . 
In addition, the refinement ( |2.12| ) leads to the expected results in simple special examples. 
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Example: Resolved Conifold 

Now, with the simple example of the resolved conifold, X = Opi (— 1) ©Opi (—1), let us 
illustrate how one can use these wall crossing formulae to obtain the (refined) D6/D2/D0 



partition functions in different chambers (|2.6|) . This discussion will be mirrored by the 



pyramid partition approach of Section 4. Before we present refined BPS invariants, let us 
briefly review the results of [0,^ for the ordinary partition functions Z{q,Q;u). Starting 
in the core region (see fig. 1) where Z{q, Q; Ccore) = 1 and applying the wall crossing 
formulae, we obtain 

n 

Z{q, Q- a+i) = - I'Qy ' (2-13) 

so that in the limit n ^ oo we recover the reduced Donaldson-Thomas partition function 

lim Z{q,Q-C^) = T\{l-q^Qy =Z'j,T{q,Q). (2.14) 



Similarly, in the chamber C^+i = [Wn+i^i 



1 

n J 5 



Z{q, Q- C.+i) = M{qf - q^QY J] (1 - q^Q'')" ■ (2-15) 

j>0 k>n 

The factor of M{q)'^ = Yl'jLii^ ~ 9"')"^"' comes from crossing the D6/D0 wall W^. In the 
limit 71 ^ oo, we now recover the full Donaldson-Thomas partition function, 

lim Z{q,Q-Cr,)=M{qfT\{l-q^Qy = ZDT{q,Q). (2.16) 

n — >oo -•- -•- 

j>0 

Now, starting with Z"^^^ {qi,q2,Q] Ccore) = 1 and applying the refined semi- primitive 
wall crossing formula ( p.l2|) we obtain the refined D6/D2/D0 partition functions in all 
chambers C^, 

n n—j+l 

^''^'(gi,?2,Q;^„+i) = n n (i-^^^r^^r^) , (2.17) 

3=1 i=l 

so that in the limit we recover a refinement of the reduced Donaldson-Thomas partition 
function: 

oo 

lim Z-f(Qi,Q2,Q;a)= n (l-^??r*9r^) ■ (2.18) 
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In the chamber = [W^^^^W^ we find 

oc oc 
i,j = l k+l>n 

(2.19) 

As in the unrefined case, the factor of M{qi,q2)'^ comes from crossing the D6/D0 waU. 
In general, when studying refined BPS invariants one encounters a family of refinements 
M5{QijQ2) = YiTj=ii^~^i ^ ^^2 ^ ^ of the MacMahon function M(q'), with different 
values of 6. For example, the refinement used in corresponds to 5 = — 1. In the classical 
limit y 1 (corresponding to qi = q2 = q) all of these refinements reduce to the ordinary 
MacMahon function M{q), while in the "opposite" limit y —>■ —1 they specialize to M{—q), 
which describes the contribution of the 0-dimensional subschemes to the DT-invariants 
of For our purposes in the present paper, it is convenient to choose a symmetric 

normalization in (|2.19 ), so that M{qi, ^2) = -^<5(?i, ^2) with 6 = 0. 

In the limit n 00 we obtain a refinement of the ordinary Donaldson-Thomas parti- 
tion function (|2.16| ), 

00 

hm (Qi, Q2, Q; = M{q,, q^f \\ [l - Q^r^^rM , (2-20) 

n — >oo \ / 

and 

^'^^'(9i,92,Q;Ci)=M(gi,g2)' n (i-Q^r^^r^) (i-Q-'^r^^r^) (2.21) 

is a refinement of the non-commutative Donaldson-Thomas partition fimction Z{q, Q; Ci). 



3. Refined = Motivic 



The unrefined wall crossing formula of Kontsevich and Soibelman [|T6[ generalizes the 
primitive (|2.8| ) and semiprimitive ( p.lO| ) cases discussed above. It encodes the degeneracies 
of BPS states in a given chamber in terms of a non-commuting product of symplectomor- 
phisms acting on the complexified charge lattice. Specifically, let 

Tr = r'^®C* (3.1) 

be an r-dimensional complex torus, where r is the rank of P, and define functions Xj 
corresponding to any 7 G P, such that X^Xy = X^+^r. Given a basis {7^} of P and 
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corresponding coordinates Xi on Tr, one endows Tr with a symplectic structure uj = 
^{li,1j)~^d log Xi A d log Xj and defines symplectomorphism^ 

: Xy ±X^)<^''^> . (3.2) 

The statement of wall crossing is that the product over all states that become aligned at 
a wall of marginal stability 

A = l[A,iu):=l[U^^^-'-\ (3.3) 

7 7 

taken in order of increasing phase of the central charge, 2(7), is the same on both sides 
of the wall. In other words, going from u = u+ on one side of the wall to w = U- on 
the other, both the BPS indices and the ordering will change but the overall product will 
remain the same: 

Yl = Yl . (3.4) 

7 7 

Geometrically, this formula arises as a specialization (namely, by taking the Euler 
characteristic) of a much more general, 'motivic' wall crossing formula |]T^ that is reminis- 
cent of the wall crossing formulae for the refined BPS invariants discussed in the previous 
section. Just as refined BPS invariants depend on an extra variable y, motivic BPS in- 
variants involve a formal variable L^/^ (where L is the motive of the affine line). Both 
refinements reduce to the ordinary BPS invariants in the corresponding limits y ^ 1 and 
LV2 ^ 1. One may therefore hope that motivic BPS invariants are precisely the refined 
BPS invariants of physics. 

There are several ways to approach such a conjecture. First, one can attempt to 
make a direct comparison of refined and motivic BPS invariants in concrete examples. 



several of which are attainable and will be discussed in |jT^. Second, the behavior of BPS 



invariants can be investigated indirectly by means of wall crossing formulae. This is our 



goal in the present paper. Using the Serre functor as in [|T6[, we pass to a "quantum" 
version of motivic BPS invariants and the wall crossing formula, and show that it agrees 
with ([2TTD . Finally, and more physically, one can try to derive the motivic wall crossing 
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See 1 16 1 and also (which we follow here in the unrefined case) for a more precise description 



of these symplectomorphisms. 
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formula via three-dimensional gauge in the presence of a graviphoton background, along 
the lines of [EB . 



The motivic DT invariants defined by Kontsevich and Soibelman are elements of 



quantum tori over a version of the Grothendieck ring of varieties [16|. The quantum 



torus in question is simply the quantization of (|3.1| ). It comprises an associative algebra 
generated by e^, 7 G F, such thalS 



^71 ^72 



,9 (71,72 > 



^71+72 



(3.5) 



and Co = 1. In particular, the generators obey the following commutation relations: 



|(7i,72> _ -|(7i,72> 1 

q ^ I 67-^-1-72 ■ 



In the classical limit, as g^/^ —1, one finds 



lim (g- 



ql/2^_l 



so that 



where 



[671 , 672 



] = (-l)^'^^'^^^(7i,72) 



^71+72 



(3.6) 



:= lim 



ql/2^_l q -1 

The Lie algebra ( |3.6| ) acts on the classical torus (pTj), and generates symplectomorphisms 
via = exp Xl^i 7^2 en7- 

Following Kontsevich and Soibelman [16|, let us introduce the "quantum dilogarithm" 
function 



E(x) := J2 



^--^(q--1)...(q--q-i 



(3.7) 



It is easy to verify that in the classical limit g^/^ —1 this function has the asymptotic 
expansion 

E(x)=exp (-^Li2(x) + ^^^ + . 

where g^/^ = —e^ and Li2(a;) = Xl^i ^ Euler dilogarithm. Moreover, the function 

E(a;) obeys the "pentagon" identity 



E(a;i)E(a;2) = E(a;2)E(a;i2)E(xi) 



(3.^ 



''' In the case of the motivic quantum torus this relation looks like 6^1672 = Lz^^i''''^) e-^_^_^_^^. 
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where X1X2 = qx2Xi and X12 = q~^^'^xiX2 = q^^'^X2Xi. This pentagon identity is the basic 
example of the motivic waU crossing formula in a simple case of two primitive charges, 71 
and 72, which obey (71,72) = 1, 

/I mot A mot /I mot a mot a mot (Q Q\ 

^1,0 ■ ^0,1 ~ ^0,1 ■ ^1,1 ■ ^1,0 • l-^-yj 

Here, := E(en7i +77172) are quantum analogs of the classical symplectomorphisms 

U^^^\ Acting by conjugation, these operators generate "symplectomorphisms" of the 
quantum torus (resp. motivic quantum torus). 

In general, the Kontsevich-Soibelman wall crossing formula for motivic DT invariants 
0™°*(7;w) says that the product of the quantum symplectomorphisms in a given sector 
does not change under wall crossing (as long as no BPS rays leave the sector): 



n^r'(^+)=n^r'(^-)- (3-10) 

7 7 

For simplicity, we assume here that each ray in this sector is generated by a single BPS 
charge 7. In this case, we have (c/. (p. 3]) above and sec. 6.2 of [|TB|) 



A^'"\u) = 1 + 2 —Al^e^ + . . . . (3.11) 

In particular, suppose 71 and 72 are primitive and consider a narrow sector containing 
2(71 +72). Then the motivic wall crossing formula ( p.lO| ) looks like 

A^MATl-y.MATM = ^r'(«-)^7°+7.(^-)^r(^-) • (3-12) 

To use this formula in practice, notice that the algebra generated by en-yi+m'y2 is filtered, 
cf. m. The lowest de gree of filtration contains the Heisenberg subalgebra 



[e^^e^J = e7,+7,, (3.13) 

with a central element e7^_|_72- Then, assuming that 0™°*(7i, w_|_) = 0™°*(7i,w_) for 
i = 1, 2 and keeping track of the coefficients of 67^+73 in ( |3.12| ), we can derive 

Al)-°H7i + 72) = (71, 72), ■ f^™°H7i)f^"°'(72) , 

where 

(71, 72)9 = 1 — zi ■ 

q2 - q 2 

This is equivalent to the refined version (|2.11|) of the primitive wall crossing formula, 
provided that we identify y ^ —q^l'^. 
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4. Refined Pyramid Partitions and the Topological Vertex 

In this final section, we interpret refined BPS invariants in terms of statistical mod- 
els of melting crystals. In the case of Donaldson-Thomas theory for non-compact toric 
Calabi-Yau 3-folds, this was done in |jT^ (unrefined) and (refined) using the topological 
vertex formalism. For other chambers of moduli space, however, rather different models 
are needed. These models are built from Calabi-Yau quivers with superpotentials; crystal 
partitions encode information about the representations of the quivers, which in turn are 



related to configurations of BPS branes. (See, for example, P,^|Tll for explicit examples 
of this relation in the Szendroi region of moduli space.) Both representations of the quiv- 
ers and BPS branes depend on a stability condition (a choice of Kahler moduli), which 
corresponds to a boundary condition for the crystal model. 

In the case of the conifold, which shall be our main example, the crystal models 
that arise are called "pyramid partitions." They essentially describe representations of 
the Klebanov-Witten quiver. By changing the boundary conditions for these crystals, 
unrefined invariants in all the chambers Cn and C„ that were described in Section 2 can 
be obtained |jT^. It turns out that by modifying the weights atoms in the crystals — 
more or less splitting weights qi and q2 across a diagonal — it is also possible to obtain 
refined invariants in all chambers. Part of our goal in this section is to describe how this 
is achieved. 

More interestingly, we find that a combinatorial transformation between partitions 



with different boundary conditions called dimer shuffling [25| corresponds very naturally 
to (refined or unrefined) wall crossing. Moreover, in the limit n oo the boundary 
conditions for pyramid partitions in chambers Cn become such that the model resolves 
into a pair of (refined or unrefined) topological vertices. This is exactly what one should 
expect after crossing the infinite number of walls to the Donaldson-Thomas chamber of 
moduli space for the resolved conifold. As the topological vertex provides a universal 
construction for the BPS partition function of toric Calabi-Yau 3-folds, we believe that 
this behavior should be quite general: a quiver-related crystal model should always resolve 
into a network of topological vertices as one approaches the DT chamber. 

We first review the pyramid partition models for the unrefined BPS invariants, and 
then proceed to refined invariants. In each case, we relate the pyramid partitions to states 
in a dimer model, in terms of which the shuffling operation (wall crossing) is most easily 
described. We then explain how both refined and unrefined models are related to the 
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topological vertex, which consists not of pyramids but of three-dimensional cubic "plane 
partitions," or three-dimensional Young diagrams. 

Unrefined invariants 




In the chamber Cm the unrefined generating function of BPS states is obtained by 
counting the melting configurations of an infinite pyramid-shaped crystal whose top row of 
atoms has length n (sometimes also called an empty room configuration, or ERC, of length 
n) @,|15[. As shown in fig. 2, this crystal has two different types of atoms, corresponding 



to the two vertices of the Klebanov-Witten quiver. The top edge of the pyramid always 
consists of n white atoms. The remainder of the pyramid is then constructed by placing 
two black atoms underneath each white atom, oriented vertically, and two white atoms 
underneath each black one, oriented horizontally.! In order for an atom to be removed 
during crystal melting, all atoms lying above it must be removed as well. The partition 
function is defined as a sum over all melting configurations [i.e. pyramid partitions) tt. 



Z{qu.,qt;Cr.)=J2l^^^''^C^^ (4-1) 



In nature, such a crystal structure, very similar to that of diamond, occurs in the minerals 
zincblende (zinc sulfide) and moissanite (silicon carbide). 
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where Wiu(7r) and Wh{'rc), respectively, are the numbers of white and black atoms removed. 
It was proven in ||2^ that this agrees with the partition function ( p.l5|) , 



Z{q^, 56; C^) - Z{q, Q- C^) = M{qf - q^QY J] (1 - ^'Q"')' , (4-2) 



k=n 



provided that one makes an n-dependent identification as in 



qw = -q'^Q 



qb = -q 



-(n-l) 



Q 



(4.3) 



Similarly, it was argued in |jT^ that to obtain the unrefined partition function in the 
chamber C„ one must sum over the melting configurations of a finite crystal configuration 
of length n, also shown in fig. 2. Then 



j=i 



if one identifies 



qw = -q'^Q 



qb 



-(n+l)g-l 



(4.4) 



(4.5) 





M 

Fig. 3: The relation between pyramid partitions and dimer states, illustrated for 
n = 2. 



odd boxes: 



even boxes: 



nr. r it 



Fig. 4: Even and odd boxes of dimers. 
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To proceed further, let us translate the above partition functions into the language 
of dimers. The partitions of a length-n pyramid correspond bijectively to the states of a 
dimer model on a square lattice with prescribed asymptotic boundary conditions. (We will 
refer to these states as partitions as well.) An intuitive way to visualize the correspondence 
(see also |^^) is to actually draw dimers on the black and white atoms, as in fig. 3. Then 
the dimer state corresponding to a given crystal automatically appears when viewing the 
crystal from above. 

As in ||2^, we have included an extra decoration on the lattices in these figures: lattice 
points are colored with alternating black and white dots. This canonical decoration carries 
no extra information, but is very useful in describing weights and wall crossing. We will 
also call squares in the dimer lattice even or odd depending on their vertex decorations. 
As shown in fig. 4, we call two dimers lying on the edges of an even (resp. odd) square an 
even (resp. odd) boT, an even (resp. odd) box with two horizontal (resp. vertial) dimers 
corresponds to a fully uncovered black (resp. white) atom in the crystal. 




Fig. 5: The weights assigned to edges of the dimer lattice of "length n," for n = 2. 
(The n = 2 ground state has been shaded in.) All vertical edges have weight 1 and 
all horizontal edges have an additional factor of {—Q)^^^^- 



One can assign weights to each edge in the dimer lattice so that the total weight of a 
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dimer partition tt, defined aM 

product of weigfits of dimerized edges in tt 

w(7T) = 

product of weiglits of dimerized edges in the ground state of the lattice ' 

agrees with the pyramid partition weight qZ'^^'^^ ■ To implement such a weighting, it 
is sufficient to ensure that the ratio of horizontal to vertical edges in every odd and even 
square, respectively, equals and q'^^ — corresponding to white atoms being removed 
and black atoms being replaced. 

Here, it is most convenient to use a weighting that is n-dependent. Vertical edges 
are always assigned weight 1. For the horizontal edges, we draw two diagonals on the 
dimer lattice, which pass through the lowermost and uppermost odd blocks in the ground 
state dimer [i.e. the lowermost and uppermost uncovered white atoms in the unmelted 
pyramid). For positive integers a, the horizontal edges 2a — 1 units above and 2a units 
below the lower diagonal are assigned weights Q)"^/^ and Q)"^^^, 

respectively, where a = means that an edge is touching the diagonal. Likewise the 
horizontal edges 2a — 1 units below and 2a units above the upper diagonal are assigned 
weights Q)"^/^ and Q)"^/^, respectively. An example is shown 

in fig. 5. For a dimer model corresponding to a length-n crystal, one can check that the 
ratios of horizontal to vertical edges in every odd block is indeed —q^Q~^ = q^, and in 
every even block the ratio is —q"'~^Q~^ = q^^ ■ Since the resulting weight function itself 
is n-dependent in terms of variables q and Q, let us call it Wn rather than w. 

We let the weight Wn be a function acting linearly on formal sums of partitions, and 
define O'^"^) to be the formal sum of all possible partitions of a dimer lattice with asymptotic 
boundary conditions corresponding to the length-n crystal. Then 



Technically, both the numerator and denominator in this definition must be "regularized." 
For a given state tt, one fixes a large box in the dimer lattice so that all dimers outside the 
box match the ground state (corresponding to an unmelted pyramid of length n); then one only 
multiplies together the weights of dimers inside this box. 
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- . i 

r 



Fig. 6: The directions in which dimers move under the shuffle S, and an example 
of shufiiing a partition of length n = 2 with odd boxes deleted. 



The operation that we claim is the combinatorial equivalent of wall crossing is de- 
scribed in as dimer shuffling. It maps partitions of length n to partitions of length 
n + 1. To define it, first consider an operation 5", which maps a dimer state tt^'^\ all of 
whose odd blocks have been deleted, to a dimer state 7?'^'^"'"^^ all of whose even blocks 
are deleted. By "deleted" we mean that any dimers forming odd (resp. even) blocks are 
removed. The operation S simply moves every non-deleted dimer one unit to the left, 
right, up, or down, according to the rules on the left side of fig. 6. We show an example 
of such a shuffling in fig. 6 as well; note that dimers carry their vertex decorations with 
them when they move. As a function from the set of {dimer partitions with odd blocks 



deleted} to the set of {dimer partitions with even blocks deleted}, S is bijective |l25|,p6 
The actual dimer shuffling operation 5" can then be defined to act on finite "subsums" in 
e(^). It maps each formal su m§ of 2"^ dimer states with a fixed set of m odd blocks (for 
any m) to the finite formal sum of all dimer states with a fixed set of even blocks in the 
obvious way: by deleting odd blocks, applying S, and filling in the missing even blocks in 
all possible combinations. Letting S act linearly on all such formal sub-sums of e(^\ it 
must, because S is bijective, send O'^"'-' precisely to G*-"^"'"^-'. 

What happens to weights under dimer shuffling? We defined our weight function 
above so that the dimers in a partition tt of a length-n model with odd boxes deleted do 
not change weight at all under the action of 5". In other words, Wn{Ti^) = Wn+i{S{Tc)). The 
only change in weights of a genuine dimer state tt under the action of S arises from the 
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We could also define shuffling, as in |25], to act on individual tt's, but this is unnecessary. 
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deletion of odd blocks and the subsequent creation of new even blocks after shuffling. An 
important lemma in (which we will refine later in this section) is that the difference 
between the number of deleted odd blocks in tt and the missing even blocks in S{tc) is 
always exactly n. Then a quick exercise shows that for a fixed tt with m deleted odd 
blocks, 

Wn(sumof ns.t.n agrees with tt) = (1 — q'^Q~^)^ ■ Wni'^') , 

w^_i_i(sumof TTS.t. TT agrees with >5(7f)) = (1 — q'^Q~^)'^~^ ■ Wn+i{S{TT)) . 

(By "agrees with," we mean aside from deleted blocks.) The ratio of these quantities is 
independent of m, immediately proving that 

^/;n(e(")) = (1 - (/"Q-^)"«;n+i(e("+i)) . (4.6) 

This is precisely the wall crossing formula between chambers C„ for the conifold. 

Formula ( [4.6|) suggests (correctly) that we can write the crystal or dimer partition 
function for a model of length n as 



w, 



(oo)^ 



Of course, the quantity Woo{Q^^^) must be the Donaldson-Thomas partition function of 
the conifold, and this relation holds because pyramid partitions of length n ^ oo effectively 
reduce to the topological vertex formalism of [p!2| , p!3[| (see also p7| , p8| ). 
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Fig. 7: The brick-like lattices around the upper and lower vertices as n ^ oo. The 
ground state state of the dimer is shaded in. As before, each horizontal edge also 
carries a weight of {—Q)^^''^ 
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To understand this relation, consider the n-dependent weighting system of fig. 5. In 
the hmit ti ^ oo, the weights of half the edges around the lower vertex (of the pyramid, or 
of the dimer model) acquire infinitely large, positive powers of q and cease to contribute 
to the partition function. Likewise for half the edges around the upper vertex. Therefore, 
the only dimer partitions around these vertices that can contribute to the length-infinity 
partition function involve dimers on edges of the brick-like lattices of fig. 7. These brick- 
like lattices, however, are equivalent to hexagonal dimer lattices, which correspond to the 
three-dimensional cubic partitions that arise in the topological vertex. 




Fig. 8: The map between the length-infinity dimer model and a pair of topological 
vertices. (The extra qi and q2 notations are for the refined case further below.) 



As argued more carefully in ||2^, any (nontrivial) configuration of the length-infinity 



dimer model can be constructed via a series of moves that amount to 1) cutting out a Young 
diagram A simultaneously from the upper and lower vertices, 2) stacking up individual 
boxes to form a cubic partition tt^ around the lower vertex, and 3) stacking up boxes to 
form a partition tt^ around the upper vertex. An example of such a dimer configuration 
and its corresponding topological vertex partitions is shown in fig. 8. By observing how 
dimers shift in these three steps and using our n oo weighting, it is not too hard to see 
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that the contributions to the partition function are (— (5)l'^lq'2 ll-^^ll 012 11-** II from step (1), 
q\'^\ I from step (2), and g'^^^tl from step (3).0 Therefore, the total partition function is 

«;oo(e(°°)) = $^ (-Q)l^i5^ilA||^+ii|A*ii^^l<j+K.-|^ 

A* ' ->> 

which is precisely the topological vertex expression for the (unreduced) partition function 
of the conifold [p!2| , p!3[| . In terms of Schur functions, the generating function for three- 
dimensional cubic partitions with a single nontrivial asymptotic boundary condition A 

is E^,?'""' = M{q)q-^\\^\\%x^{q-P) = M(q) ^"^11^11%;,* (q1/2^ ^3/2^ ^5/2^ Thus, as 
expected, 

w^Q^^^) = Z{q,Q-C^) 

= M{qfY.^-Qfs^{q-P)sA<l-n 

00 

= M(q)2 J] (1 - Q^-l/2gfc-l/2Q) (4.7) 

00 

^M{qfl[{l-q^Qy. 

Refined invariants 

We now come to the crystal melting models for refined invariants. For the conifold, 
we describe models to compute the refined partition functions in all chambers Cn and 
Cn- We will first present the formulae in terms of melting crystals, and then prove them 
while discussing their relation to dimer shuffling, refined wall crossing, and the the refined 
topological vertex. 

At the level of crystal models, one must draw a series of diagonals on the pyramid 
partition, and interpolate weights between the variable gi on one side of the diagonals and 
q2 on the other. To be more specific, consider the pyramid of length n = 1, corresponding 
to the Szendroi chamber Ci. On this crystal model, we draw a single diagonal as shown 
in fig. 9; we assign white atoms above the diagonal a weight g^, white atoms below the 
diagonal a weight g~, and white atoms on the diagonal itself a weight (q^Qw)^^"^- black 
atoms are assigned weight g^. Letting zi;+(7r), w~{7t), and w^(7r) be the numbers of white 

We use conventional notation for Young diagrams and three-dimensional cubic partitions; A* 
is the transpose of the diagram A, the rows of A have lengths Xi, |A| = '^Xi is the number of 
boxes in A, ||A|p = ^ A^ , and |7r| is the number of boxes in a three-dimensional partition. 
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Fig. 9: Weights of atoms for the refined partition function in chamber Ci. 



atoms above, below, and on the diagonal, respectively, in the partition tt, and identifying 
(it = -QiQ~^, Qw = -?2Q"^ and Qb = -Q, we find 

71" 

= Z^^^iqi, Q2, Q; Ci) (^4_g^ 



To generalize to the length-n pyramid, we draw n diagonals, as in the left half of 
fig. 10. It is more natural to work directly in terms of the variables qi, q2, and Q. We 
assign weights —qiQ~^ (resp. —qi ^^Q) to the white (resp. black) atoms above all the 
diagonals and weights ~q2Q~^ (resp. —q2 ^^Q) to the white (resp. black) atoms below 
all the diagonals. The diagonals themselves intersect white atoms; we assign the same 
weight to all the white atoms on a single diagonal, interpolating between —q^ ^q^Q 
on the uppermost diagonal and — q'2 ^ on the lowermost (multiplying by qiq2 in 
each intermediate step). Similarly, black atoms lie between diagonals, and we assign them 
weights ranging from —q^ ^q2^Q directly below the upper diagonal to —qi^ ^2 
directly above the lower diagonal. Multiplying together the weights of all atoms removed 
in a given partition tt and summing these quantities over partitions, we obtain the expected 

00 

z-f(Qi,Q2,Q;a)=M(Qi,Q2)' n(i-5r^9r^Q) n (i-^r'^r^Q"')- (4.9) 

i,j = l i>i,i>i 
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Fig. 10: Refined weights of atoms for chambers Cn and Cn+i, with n = 3. 

For chambers Cn+i, the finite pyramid of length n can also be split by n diagonals, 
as shown in the right half of fig. 10. If one assigns weights such that (1) when qi ^ q 
and q2 q white atoms have weight —q^Q and black atoms have weig 
(2) when moving up one step, either on or inbetween diagonals, the absolute value of the 
power of (j2 (resp. qi) decreases (resp. increases) by 1; and (3) the assignment is symmetric 
about the middle diagonal(s) of the crystal, the resulting partition function is precisely 

^'^^'(9i,?2,Q;c'„+i)=M(Qi,Q2)' n (i-^r^^r^Q)- 

i>l,j>l 
i + j <™ + l 

For the remainder of the section we return to the infinite pyramid of length n, gen- 
eralizing the previous unrefined discussion to refine the connection between shuffling, wall 
crossing, and the refined topological vertex (and to prove formula (|4.9| )). We first observe 
that in order to equate refined pyramid partitions and their weights with states (configu- 
rations) of a dimer lattice, we can use almost the same n-dependent weighting described 
in fig. 5. Now, for positive integers a, the horizontal edges 2a — 1 units above and 2a units 
below the lower diagonal are assigned weights ^''^^(— Q)"^/^ and q2 ^''^^(— <5)~^^^, 
respectively. Likewise the horizontal edges 2a — 1 units below and 2a units above the 
upper diagonal are assigned weights q^"" ^^^"^ {—Q)~'^^'^ and q^ ^■*^^(— Q)""*^/^. See the 
example in fig. 11. 
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Fig. 11: Refined weighting of the length-n dimer, for n = 2. (The color coding 
here differs from that in fig. 5, to emphasize the difference between q\ and 52 •) 

As in the unrefined case, the weights of dimers which are not part of deleted odd or even 
blocks do not change during dimer shuffling S, due to our n-dependent weighting. In order 
to understand the behavior of the deleted blocks, we observe that the shuffling S removes 
exactly one (deleted) odd block from each of the n diagonals of a dimer configuration of 
length n. Moreover, the remaining (deleted) odd blocks are mapped to deleted even blocks 
with exactly the same weights — if (for instance) they were above all the diagonals, then 
they remain above all the diagonals. These statements can be proved with careful counting 
arguments, considering the number of dimers on and around each diagonal in an arbitrary 
configuration before and after shuffling. The result is that when the actual shuffling S 
maps a formal sum of states tt agreeing with a fixed odd-deleted state tc on all but their 
odd blocks to a formal sum of states agreeing on all but their even blocks, the weight of 
this formal sum changes by exactly ni+j=n+i(-'- ~ 5i ^^2 therefore, 

«^n(e(-)) = n (1 - ^r^^r^Q"') ■ «^n+i(e(-+i)) . 

i > 1 , J > 1 
i+j = n + l 

This, of course, is the refined wall crossing formula for chambers Cn- 
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Fig. 12: Neighborhoods of the refined upper and lower vertices as n ^ oo. 



The crystal-melting or dimer partition function of length n can now be written as 

oo 



The last term, Woo(0*-°°-*), is obtained from a slightly modified version of the refined 
topological vertex of 0. To see this, observe that as n ^ cxd the neighborhoods of the 
upper and lower vertices of the dimer lattice still reduce to effective brick-like lattices, 
now shown in fig. 12. In terms of three-dimensional cubic partitions, states of the length- 
infinity dimer are again created by 1) cutting out a Young diagram A simultaneously from 
the upper and lower vertices, 2) stacking up individual boxes to form a cubic partition 
Tv^ around the lower vertex, and 3) stacking up boxes to form a partition around the 
upper vertex. The creation of the Young diagram A comes with a fairly simple weight 



(— (5)'^'q'i Qi^^^ " • However, both the upper and lower "room corners" are now split 
along a diagonal, as shown in fig. 8. In the case of the lower corner, boxes stacked below 
the diagonal come with weight q2, those above the diagonal with weight gi, and those that 
the diagonal intersects have weight {qiq2)^ ■ The situation is reversed for the upper vertex. 
The generating function for such three-dimensional cubic partitions with one asymptotic 
boundary condition A is (for example, at the lower vertex) 



(qi) 



Q2 



(9192)- 



M{qi,q2)q2 



(4.10) 
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with M(gi, = Hi j=i(1^9i ^ I2 ^) ■ Therefore, the length-infinity pyramid partition 
function is 



= M(gi,g2)'$^(-Q)l"I^A(?2"')«A*(9r) 



A 

00 



Note that expression (|4.10| ) differs only slightly from the refined topological vertex used 
in (with the boundary condition A placed along an "unpreferred" direction) . The differ- 
ence comes from our symmetric choice of normalization, as discussed in section 2. In the 
diagonal is assigned to q2 rather than {qiq2)^ , resulting in the fact that the refined MacMa- 
hon function appearing in the analogue of ( |4.1U| ) is not M{qi, ^2) = Y[Tj=i(^~'ii ^ ?2 ^ )~^5 
but rather "~ Qi~^Q2)~^- '^^^ refined A- mo del (Gromov-Witten) partition functions 
calculated with the refined vertex are always normalized by the prefactor M{qi,q2)^, so 
in many previous calculations this has made no difference. 
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